arXiv:nucl-th/0112082vl 30 Dec 2001 


Nuclear magnetization distribution and hyperfine splitting in Bi S2+ 

ion 

R.A. Sen’kov and V.F. Dmitriev, 

Budker Institute of Nuclear Physics, Novosibirsk-90, 630090, Russia 

February 7, 2008 


Abstract 

Hyperfine splitting in Bi 82+ and Pb 81+ ions was calculated using continuum RPA approach with 
effective residual forces. To fix the parameters of the theory the nuclear magnetic dipole moments of 
two one-particle and two one-hole nuclei around 208 Pb were calculated using the same approach. The 
contribution from velocity dependent two-body spin-orbit residual interaction was calculated explicitly. 
Additionally, the octupole moment of 209 Bi and the hfs in muonic bismuth atom were calculated as well 
in the same approach. All the calculated observables, except the electronic hfs in 209 Bi, are in good 
agreement with the data. We argue for more accurate measurement of the octupole moment and the 
muonic hfs for 209 Bi. 


1 Introduction 

High experimental precision attained in the laser spectroscopic measurement of the ground-state hfs in 
hydrogen-like 209 Bi 82+ JlJ stimulates considerable theoretical activity in this field (see, e.g., |n] |nj and 
references therein). 

A E e h x / s er - = 5.0840(8) eV. (1) 

The first quantum mechanical calculation of the hfs was made by Fermi in 1930 ||]. He treated the hfs 
with a non-relativistic formalism in studies of alkali atoms, where he derived approximate generalizations 
from the hydrogenic case. In the case of 209 Bi 82+ the formulae obtained by Fermi give a value about 2.75 eV, 
which is 46% smaller than the experimental value (fl). Relativity can be taken into account by multiplying 
non-relativistic value with a relativistic correction factor A(aZ) ||. In bismuth case Z = 83 and factor 
A(aZ) = 2.125. The relativistic generalization of the Fermi formulae holds for point-like nucleus and gives a 
hfs for 209 Bi 82+ of about 5.84 eV, which is 15% larger than the experimental value ([!]). 

The effect of spatial distribution of the nuclear charge was analyzed in pioneering works by Rosenthal, 
Breit Q] and others |^],|(| and it is sometimes called ”Breit-Rosenthal effect”. The correction factor for the 

nuclear charge distribution can be written as (1 — <5), where S is a small number which depends mainly on the 

I/2 

root-mean-square (rms) radius of the nuclear charge distribution, (r c 2 ) . By assuming a uniform spherical 

symmetric charge distribution and using the experimental value (rc 2 ) 1 ^ 2 = 5.519(4) fm JtJ] for the nuclear 
charge distribution, S can be calculated to be 0.0110 in the case of 209 Bi 82+ Q. The relativistic hfs energy 
splitting in 209 Bi 82+ for a uniformly charged nucleus becomes than about 5.20 eV, which is 2% larger than 
the experimental value (|l|). 

There are two principal corrections to the Fermi-Breit formulae that is necessary to take into account: 
the magnetic moment distribution within the nucleus i m and radiative corrections JTi],|h] [h]]. Both of 
these corrections are of comparable magnitude. 

The correction for an extended nuclear magnetization was first studied in an innovative work by A. Bohr 
and V. Weisskopf H[nJ and is also called ”Bohr-Weisskopf effect”. Recently, three approaches for theoretical 
determination of the Bohr-Weisskopf effect have been used. The simplest of these approaches is based on a 
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solution of the Schrodinger equation for a nucleon in a mean-field potential, the solution gives the distribution 
for the unpaired nucleon in the nucleus and the distribution is used to determine hyperfine splitting. This 
approach has been used by Shabaev et al. 0] and by Gustavsson et al. 0]. A slightly more sophisticated 
approach, giving equivalent results for the hfs in leading order, is the ’’dynamical proton model” (DPM), 
where the odd proton of the bismuth nucleus is treated as a Dirac particle bound in a mean-field potential. 
The first order hfs in hydrogen-like Bi is then given by a photon exchange between the electron and the 
proton. DPM was introduced by Labzowsky et al. [0 1. The third and more complete approach is a rnany- 
body calculation with the use of ’’dynamical correlation model” (DCM). This approach is the only calculation 
of the Bohr-Weisskopf effect which includes many-body contributions and Tomaselli et al. have used it for 
studies of several systems 0 - 0 - 

In addition to the dominant electrostatic and hyperfine interactions with nucleus, the electron also inter¬ 
acts with radiation field, an interaction described by QED. The leading QED corrections originate from the 
one-loop self-energy and vacuum-polarization effects. The one-loop QED effects for hfs have been calculated 
by different groups and the results are consistent 0,0 0. 

In this work we consider a contribution of the magnetic moment distribution within the nucleus to the 
hydrogen-like ion hfs. Our approach is close to the one used in 0]. This microscopic approach is based 
mainly on the Migdal’s theory of finite Fermi system (FFST) 0] (see also 0 0). Essentially, this approach 
is equivalent to calculation of the core polarization effects with the use of an effective interaction. The first 
calculation of the core polarization goes back to fifties 0] and this early development was summarized in 0]. 
Modern calculations, however, differ significantly in the size of the single-particle space included into sum 
over intermediate particle-hole states. For spherical nuclei all single-particle space, including continuum is 
included into calculation. This is one difference between our approach and the DCM where one has to use 
the restricted single-particle space. Another difference between our work and 00 is that RPA accounts 
for ground state correlations represented graphically by backward going loops. And, finally, we use explicitly 
the correction to nuclear electromagnetic current due to velocity dependent interaction. This correction was 
absent in 00 since they use velocity independent forces. Previously, this correction was discussed in p(i| 
in connection to the Bohr-Weisskopf effect in bismuth muonic atom. It was shown that this correction is not 
negligible. Here we extend this analysis to the hydrogen-like ions around 208 Pb nucleus. 

2 Basic Equations 

In FFST the effect of the core polarization is described by introducing an effective single-particle operator 
(vertex) M satisfying the equation 0] 


M M = e e ff M 0 + F A(w) M(u), (2) 

where M$ is in our case the bare single-particle MJ operator, A(w) the polarization operator of a particle-hole 
pair, 

A iX 2^3 ^4 = y V3 (e + ui) G „ 2 „ 4 (e - -w), (3) 

where G uu ' (e) is the single-particle Green function. F is the amplitude of quasi-particle interaction. Equation 
(]2j) is written in the re-normalized form, i.e. after elimination of a regular part in the Green function G^'le). 
The regular part of the G-function corresponding to admixture of three- or more particle states produces 
the effective charges e e //, and re-normalizes the interaction amplitude F in (0 Following [0, we choose 
e e ff as a constant independent of the particular state \v) 1 but different for the spin and orbital parts of the 
operator M 0 : 


(1 - C) (M ,% + C (M§)n , 

for protons, 

(4) 

(1 - C,) (M 0 S )„ + C (M 0 S ) P , 

for neutrons, 

M i 1 (1 - CO (M l 0 ) p , 

eff M 0 _ | C (Mq) p , 

for protons, 
for neutrons. 

(5) 
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The constants £ in (||) and (|B|) were taken as equal for both protons and neutrons since the main deviation 
of the magnetic moments from single-particle values is in the isovector part of the MJ operator. 

The effective interaction is F = F ss + Fi s , where 

F ss = C (g + g' Ti • t 2 ) cr 1 • <x 2 <5(n - r 2 ) (6) 

is the Migdal-type spin-spin zero-range interaction and 

Fu = Crl (k + /c'n ■ r 2 ) ■ (Pi - p 2 ) • (o^ + <x 2 ) x Vi<5(ri - r 2 ) (7) 


is the spin-orbit one. 

The M J vertex M jm can be written in coordinate representation in the form 

3 

= vi ( r ) T JM > 

i—1 

where we have introduced the complete set of linear independent tensor operators 

Fjm = Tjm = a 

Fjm = H L • Y »(n) + Y^( n) ■ L) , 

where Yj M (n) being the vector spherical harmonic pfS| |. 

The values of the MJ moments can be expressed in terms of matrix elements of the vertex (|j): 

IJ-J = (vm = I\ M J0 {u) = 0) \um = I) = ( ^ nlI \ u ’( r ) M-0 


( 8 ) 


(9) 


( 10 ) 


where v = ( nllm ) is the set of nucleon quantum numbers, and = ^Z 2 7 2 ||Tj I ' ) ||Zi/i^ is the reduced 

matrix element of the tensor operators (|)J). 

Reducing the angular and spin variables in equation ([|) we obtain in coordinate representation the system 
of integrals equations 

3 r oo 

w l {r) = e e ffWo l (r) + V"' / dr'O 1 -* (r, r'\u>) vP (r'), (11) 

j = l J 0 

where we have introduce w l (r) = rv l (r). The kernel of the integrals equations (0 is the sum of two terms 

e y = + eg. (12) 

For the spin-spin interaction (j(S]) 

= Cff rAy(r,r'\us)r', (13) 

where g = g + g 1 t± ■ t 2 and the sun over isospin variables is assumed. The polarization operator A’j (r, r'\ui) 
can be calculated in terms of the Green function Gij(r,r' |e) of the radial Schrodinger equation: 

A j{r,r'\u)) = 2 j +l Y ( r ') 

x [Ghj 2 (r,r'\e vi - w) + G hh (r,r'\e Vl + u;)], (14) 

where k„ is the occupation number of the level v = ( nljm ), and R v {r) is the single-particle radial wave 
function. The Green function can be calculated using two linear independent solution of the radial Schrodinger 
equation. This method allows one to use all of the single-particle spectrum |2^], [Q. For the spin-orbit 
interaction ([?]) expression in terms of A l j is more complicated (see p0[). 
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3 Hyperfine Splitting and Static MJ-moments 

The ground-state hyperfine splitting of hydrogen-like ions is conveniently written in the form 0, 0 
AE hfs = ^a(aZ ) 3 —— | 1 m e c 2 A{aZ){ 1 - £)(1 - e) + A Eqed- 

6 hn Trip 21 


( 15 ) 


Here a is the fine-structure constant, Z is the nuclear charge, m e is the electron mass, m p is the proton mass, 
/r is the nuclear magnetic moment, /ijv is the nuclear magneton, and / is the nuclear spin. A(aZ) denotes 
the relativistic factor [[| 

A{aZ) = . 1 — -r, where 7 = y/l - (aZ) 2 . (16) 

7(2y - 1) 

<5 is the nuclear charge distribution correction, e is the nuclear magnetization distribution correction (the 
Bohr-Weisskopf correction) |J, and A Eqed is the QED correction. To obtain the first part of Eq. ©, one 
can start from the magnetic interaction of an electron and a nucleus 


H int = —V- / d 3 rd 3 


„, j(r)-J(r') 
|r —r'l 


where j(r) and J(r') are the electromagnetic current densities. The hfs can be written as 

AEhfs = AEi + A Eqed, 

where AE± is the diagonal matrix element of 0-th spherical component of the vector vertex Ti M (r), 

AEi = ( I'm = I\ViQ\vm = I). 

The vertex T 1/4 is the solution of equation (j|) with the bare vertex T ] 1 derived from (|l7|), 

AeK 2/ + 1 


(17) 

(18) 

(19) 


F° = 

lfi j + 1 2 1 


Miflnij - Y [&>( a ) T il\ a ) + 9i(a ) T^(a) ft %j(r a ) 

l a 

[^(«) T 1^(«) - \j\ 9s{a)T[^{a) fi^ra) j . 


( 20 ) 


In this expression the sum is taken over the nucleons, K = (l — j) (2 j + 1) for the atomic states with j = l ± |, 


n<°> = 

rdj 


dr' fnij{r')g n ij{r'), 


= f dr' frdj {r')g n ij (r'), 
Jo 

Q nij( r ) = 4 [ dr'r' 3 f n ij {r')g n ij (r 1 ), 
r J o 


( 21 ) 

( 22 ) 

(23) 


where f n ij(r ) and g n ij{r) are the upper and the lower components of the Dirac radial wave function. The 
vertex is the bare operator of the nuclear magnetic dipole moment. The bare operator of the nuclear 
magnetic multipole moments has a standard form f27]] 


M° jm = V^l^Y rJ a 1 (^( a ) T «( a ) +9i{a) j-^Y T JM( a )V 


(24) 


The vertex DA is Ml part of the interaction ( 0 ) averaged over atomic state with the quantum numbers 
n,l,j. Higher moments of the interaction ( 0 ) do not contribute to hfs for the electron ground state with 
3 = 1 / 2 - 
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Table 1: Interaction parameters 


C (MeV- fm 3 ) 

r 0 (fm) 

K 

k’ 

9 

9' 

300 

1.16 

0.21 

-0.12 

0.633 

1.0133 


The same equations (2j) can be used to calculate nuclear magnetic multipole moments At jm with the 
bare vertex (p4|). The Eqs. (||^), (|d[) obtained using the standard electromagnetic current density for free 
nucleons. However, due to velocity dependence of the spin-orbit interaction (R), there are the corrections to 
electromagnetic current density. The corrections were discussed earlier |po|,|27|. They produce an additional 
contributions to the bare vertexes Ad ° JM and T° and can be presented in the similar form as (p(i|) , (24). The 
corresponding equations can be found in the Appendix A. 

It is worth to note, that the expression for hfs (^0|) is a sum of three terms. The first, being proportional 
jj, corresponds to a point-like magnetic moment distribution. It does not contribute to the correction 
e. This term gives the main contribution to hfs being larger than two other terms by about two orders of 
magnitude. The effects of nuclear structure are not essential here, they are hidden in the value of the nuclear 
magnetic moment which is well known from the experiment. Two other terms are, just, the correction for an 
extended nuclear magnetization. But, only the second term is proportional to the magnetic moment density, 
while the third one is not. 

The kernel ([E|) of the Eq. (jTlj) was calculated using partially self-consistent mean field potential [^8) . The 
potential includes four terms. The isoscalar term is the standard Woods-Saxon potential 


Uo(r) 


V 


1 + exp 


r-R ’ 


(25) 


with the parameters V = 52.03 MeV, R = 1.2709A 1//3 fm, and a = 0.742 fm (28). Two other terms Ui s (r), 
and U T (r ) were calculated self-consistently using two-body interaction Eq.(Q) for the spin- orbit part of the 
potential, and the Migdal-type isovector interaction 

F t = f'C(r i • r 2 )<5(ri - r 2 ) (26) 

for the isovector part of the potential, with f = 1.075 (2^]. The last term is the Coulomb interaction that 
was calculated for uniform charge sphere distribution with Rc = 1.18A 1 / 3 . This potential produces for 209 Bi 

■j- 1 2 

nucleus a charge density with rms. (r c 2 ) = 5.52 fm, which is in fair agreement with the measured value 

5.519(4) fm §. 

The parameters of the interactions Eqsdjoj.jT]) used in calculations are listed in Table 1. This is a standard 
set of the parameters used in lead region plf. 


4 Effects of the core polarization 

Qualitatively, the effects of the core polarization are well understood. The spin-spin interaction is repulsive, 
therefore, it creates negative core response decreasing the spin contribution to the magnetic moment. Fig. 
1 shows the core polarization effects on the spin part of the magnetic moment in 209 Bi. The minimum near 
nuclear surface is due to transitions to the partially filled upper level lhg/ 2 of the spin-orbit doublet. The 
transitions to higher states, including continuum, produce a uniform decrease of g s near r = 0. The core 
polarization effects are more pronounced in the octupole moment of 209 Bi. Fig. 2 show how the original r 2 
dependence is modified by the transitions over Fermi surface. The number of the transitions is greater than 
in case of magnetic dipole moment. Here the transitions with AN = 2 are involved as well, therefore, the 
relative effect of the polarization is greater than for the magnetic moment. 

Although the core polarization increases absolute values both of the magnetic moment and the octupole 
moment, this increase is not enough to explain the data. Table 2. shows the contributions to the magnetic 


5 












r [ fm] 


Figure 1: Effects of the core polarization on the 
spin part of the magnetic moment in the 209 Bi. 
Dashed line shows unrenormalized g s = 2.79. Full 
line shows renormalized g s as a function of r. 



r [ fm] 


Figure 2: Effects of the core polarization on the spin 
part of the octupole moment in 209 Bi. Dashed line 
shows unrenormalized Q s = r 1 2 dependence. Full line 
shows renormalized f l s as a function of r. 



M v 

A M u 

M 

AM 

M + AM 

Exp. 

dipole 

2.62 

0.49 

3.28 

0.30 

3.58 

4.110(4) 1 

octupole 

-14 

-9 

-44 

-4 

-48 

-55(3) 2 


Table 2: Contributions to the magnetic dipole and octupole mo¬ 
ments in the 209 Bi, 


dipole and octupole moments in 209 Bi from M° and AM 0 and the effects of the core polarization when 
Qs = Q = 0. The remaining contribution comes from higher orders configuration mixing, and from other 
possible velocity dependent forces not included into consideration. This contribution is expected to be small. 
It can be taken into account phenomenologically, via the effective charges ( s and Q. Fitting the value of the 
209 Bi magnetic moment only is not enough to find both £ s and Q. To obtain the region of the allowed values 
for and Q we calculated the magnetic moments for all four nuclei lying near doubly-magic nucleus 208 Pb. 
They are 20 ' Tl, 207 Pb, 209 Bi, and 209 Pb. Our suggestion is that for these four nuclei the difference between 
their effective charges is small and can be neglected. A calculated magnetic moment is a linear function of the 
giromagnetic ratios g s ,gi- So, on the ( s , Q plane, all points where the calculated magnetic moment is equal to 
its measured value, lie on a straight line. In Fig. 3 we plotted the lines corresponding to fits of the magnetic 
moments for the set of nuclei, mentioned above. The data were taken from Ref. ]3l]]. The accuracy of the 
measurement is high, and two lines corresponding to the experimental values g ± Ag almost coincide within 
the scale of the figure. In an ideal theory all four lines should cross in the same point. This is not the case 
for our model. The allowed values for £ a and Q form the whole region which is shown in Fig.3 by the shaded 
triangle. The “best” values of £ s and Q were found by minimization of a mean squared deviation of the 
calculated magnetic dipole moments for these four nuclei from the data. They are Q = —0.074, ( s = 0.030. 
These values are small enough in line with the above suggestion. The size of the shaded triangle in Fig. 3 
can be used to estimate the accuracy of our approach. In our model the uncertainties in calculation of the 
hfs were found by varying ( s and Q within this shadow triangle. 

The resulting distribution of the magnetic moment density in 209 Bi is shown in Fig. 4. The magnetization 
density is peaked near nuclear surface both for the bare contribution of the unpaired proton, shown by a 
dashed line, and for the total density shown by full line. The peak position in the total magnetization density 
is shifted slightly to smaller r. This is related directly to the decrease of g s inside the nucleus (see Fig. 1). 
The decrease leads to enhancement of the magnetic moment in this region due to opposite sign of the spin 

1 The value was taken from Ref. ( 2£ ] 

2 The value was taken from Ref. [ 3C ] 
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Figure 3: Region of allowed values for Q and £s with core polarization effects. 

Table 3: Theoretical values of the ground-state hfs and the experimental results. The major contribution to 
uncertainties of the theoretical values originates from the uncertainties in the Bohr-Weisskopf effect. 


Ehfs in eV 

209 Bi 82+ 

20 7 p])81 + 

Eq 

A Ebw 
A-Eqed 

j-ttheory 
hf s 

rpexper. 
hf s 

5.191(5) 

-0.050(- 2 3 o ) 

-0.030 

5.111(- 2 3 o)(5) 

5.0840(8) 3 

1.274(3) 

-0.045(+ 4 ) 

-0.007 

1.222(+ 4 )(3) 

1.2159(2)0 



e for 209 Bi 82 + 

e for 207 Pb 81+ 

Our work 

Shabaev et al. |ll|] 
Gustavsson et al. Q 
Tomaselli et al. fl5|| 
Experiment 

0.0095(+3 7 8 ) 

0.0118 

0.0131(26) 

0.0210(17) 

0.0147(11) 

0.0353(- 33 4 ) 

0.0419 

0.0429(86) 

0.0289(15) 

0.0397(25) 


part of the magnetic moment relative to its orbital contribution for /i g / 2 level. The shift of the peak results 
in some decrease of the magnetization rms-radius. The calculated magnetization rms-radius is practically 
insensitive to the values of Q and £ s within the above range. It is equal to (r^) 1 / 2 = 5.86 fm. Fig. 5 gives 
similar figure for the octupole moment density distribution in 209 Bi. In Figs. 4,5 one can see that the core 
polarization effects for the octupole moment are really larger than for the magnetic moment. 

The results obtained for hfs are listed in the left part of Table 3. Here Eo corresponds to the relativistic 
Fermi formulae plus charge distribution correction. The uncertainty in the value of Eq originates mainly 
from the uncertainty of the measured magnetic moment: ^ = 4.110(4)^ for 209 Bi and /.t = 0.5925(6 )/j,n 
for 207 Pb ^9|. The uncertainty of the charge distribution correction is small and can be neglected. The 
uncertainties in Bohr-Weisskopf effect A Ebw come from the uncertainties in the effective charges £ s and Q. 
They were determined using the shadow triangle shown in Fig. 3. They are not symmetric since the “best” 
values of and ( s are lying close to 209 Bi line in Fig. 3. In order to compare our calculations with the data 
we need A Eqed- It was calculated in several papers and the results are consistent |lI],|TGj [l^]. In |l^] it was 
shown that A Eqed is insensitive to details of the magnetization distribution. They give for A Eqed the 
value A Eqed = —0.0298(3) eV for 209 Bi 82+ and A Eqed = —0.00726(7) eV for 207 Pb 81+ . In our calculated 
values of the ground-state hfs E the uncertainty in the first brackets originates from A Ebw and the 
second one originates from A Eqed ■ The right part of the Table 3 gives our results for Bohr-Weisskopf effect 

3 The value was taken from Ref. 

4 The value was taken from Ref. 
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r [ fm] 


Figure 4: Magnetic moment density distribution in 
209 Bi. A contribution from a bare unpaired proton 
is shown by a dashed line. Full line shows the total 
magnetic moment density including core polariza¬ 
tion. 



Figure 5: Octupole moment density distribution in 
209 Bi. A contribution from a bare unpaired proton 
is shown by a dashed line. Full line shows the total 
octupole moment density including core polariza¬ 
tion. 


in comparison with some other theoretical calculations. While our result for the ground-state hfs of 20, pb 81 + 
is in good agreement with the data, the result for 209 Bi 82+ lies out of the data. The calculated e at its upper 
limit e = 0.0103 differs by four standard deviations from the measured value e = 0.0147 ± 0.0011. Actually, 
e is not directly measured value. It was obtained from the measured E < f^ p s er ' using the calculated values of 
Eq and AEqed- The uncertainty in e comes mainly from the uncertainties in these calculated E 0 . In other 
calculations only in [|I5| the particle-phonon coupling was accounted in scope of the DCM. In |ll|, and [Q, 
and others (see Refs, in Jl2[) the unpaired particle was treated as an independent particle moving in a mean 
field potential. In this approach, the difference in e reflects more the sensitivity to the choice of a particular 
mean field potential. It is interesting to note that e calculated in jl5) lies almost on the same distance from the 
experimental value, but on the other side compared to our value. This difference may be attributed to ground 
state correlations that contribute significantly to transition probabilities at small excitation energies |34). In 
order to check whether our results is sensitive to a particular parameterization of the effective charges we 
made calculation of the hfs in 209 Bi 82+ treating all the giromagnetic ratios gi(p), gi(n), g s (p), and g s (n) as free 
parameters. Their values were obtained by fitting the magnetic moments of four nuclei mentioned above. The 
magnetic moments were calculated together with the core polarization effects and the /s-corrections. With 
these new parameters we calculated hfs for 209 Bi 82 + and obtained exactly the same result A Ebw = —0.050 
eV. This situation seems to be rather general. We changed different parameters entering in our theory, 
including g and g' , although g' is fixed by the position of Gamow-Teller resonances |3^]. But, as soon as the 
magnetic moment of 209 Bi is fitted, the correction A Ebw becomes close to the value cited in Table 3. 

In addition to the magnetic moment, for 209 Bi there are other observables related to the magnetic prop¬ 
erties. In Table 4 we summarized our results including there our calculations of the octupole moment of 
209 Bi and the hfs in muonic 209 Bi. To demonstrate the relative importance of the core polarization effects 
for different observables we made additional calculation switching off the core polarization and keeping the 
same parameterization of the effective charges via C,i and £ s . One can see that although the effects of the core 
polarization are not very significant for the electronic hfs, for other observables, like the octupole moment, 
they are rather large. All the observables except the electronic hfs in 209 Bi are in good agreement with the 
data. However, the experimental uncertainties in case of the octupole moment and the muonic hfs are rather 
large. It would be very desirable to reduce them in order to see whether the discrepancy in the electronic hfs 
in 209 Bi were pronounced in these observables as well at higher accuracy of the data. 

In summary, using continuum RPA with effective residual forces we calculated the distribution of the 
magnetic and the octupole densities for 2Q9 Bi. Additional contribution from velocity dependent spin-orbit two- 
particle interaction was calculated explicitly. The parameters of the theory were fixed by fitting the magnetic 
moments of four one-particle and one-hole nuclei around 208 Pb. Basing on these results we calculated the 
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Table 4: Relative importance of the core polarization effects for different observables. 



Without core p. effects 

With core p. effects 

Experiment 

A E bw for 209 Bi, eV 

-0.046(- 3 3 9 ) 

-0.050(- 2 3 o ) 


E hfs for 209 Bi, eV 

5.115(- 3 3 9 )(5) 

5.111(- 3 0 )(5) 

5.0840(8) 

A Ebw for 207 Pb, eV 

-0.035(+ 3 7 ) 

-0.045 (+ 2 4 i) 


E hfs for 207 Pb, eV 

1.232(+ s 3 7 )(3) 

1.222(+ 4 )(3) 

1.2159(2) 

A Ebw for muonic 209 Bi, KeV 

-l-73(+iL) 

-1.86(-^) 


Ehfs for muonic 209 Bi, KeV 

4.76( + -® 9 )(6) 

4-63(; 4 °)(6) 

4.44(15) | 

Octupole moment of 209 Bi, pN-im 2 

-42(+3o) 

-55(4) 

-55(3) 


hfs in hydrogen-like 209 Bi 82+ and 20 'Pb 81+ ions, the octupole magnetic moment in 209 Bi, and hfs in muonic 
atom of 209 Bi. Except the electronic hfs in hydrogen-like 209 Bi 82+ all other results are in good agreement 
with experiment. 
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A Spin-orbit corrections to the bare vertices 

The vertex Mjm is the solution of equation (j§|) with the bare vertex (|24|). The Eqs. (20), ( |24| ) obtained 
using the standard electromagnetic current density for free nucleons. However, due to velocity dependence of 
the spin-orbit interaction (7), there are the corrections to electromagnetic current density. The corrections 
were discussed earlier [|o[|2^]. They produce an additional contributions to the bare vertexes M° jm and T^, 


AM 0 f M = V^l m ^Cj2 


-(K + n')r J a 1 p p (r a ) + (k - k') 


1 jdp n (r a ) 


2J + 1 


dr 


Tjm ( a ) 


+ (ft — k ') 

AM% = V4^Jp N =^Cj2 


(2J + 1) V J + l 


J jdpn{r a ) T ( 2 ) / \ 

' n : 1 JM\ a h 


dr 


-{K-n')r J a 1 p p (r a ) - (k- k') 


1 j dp p {ra) 


2J + 1 a dr 


Tjm ( a ) 


— (k — k')- 


J r j dp P (r a ) (2) , , 

■ T " dr 


Ar 0p= 4e/\ 2J + 1 
j + l 2 1 


(27) 


(28) 


(2 J + 1) V <7 + 1 

A ^lp nlj ~ - (V + K ')Pp( r a) ~ - K ') r a dPn ^ a ^ ^ (r a ) 




rfiV) + 


l( K _ K>a ^n«(r a ) 


(K + K')p p (r a ) + - K’)r a df>n ^ \ ft iy(r a ) 




(29) 


5 The value was taken from Ref. nr 
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Ar?£ = AeK 21 + 1 


3 + 1 2/ 


A Ad?” Q nlj - V^VN = ^T L C'Y^ - f(K - K')p p (r a ) + ^( K ~ K ') r a dpp d ^ a ' > ) Q nlj( r a) 

n. *- \ ' 




dr 


ifi’W + 




dr 


+ f-(« - K')p p (r a ) + - n')ra dPp ^ r a " > \ fiiy(r a ) 


Ti ( ?(o) 


(30) 
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